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ABSTRACT 

Massive central objects affect both the structure and evolution of galactic nuclei. Adiabatic growth of 
black holes generates power-law central density profiles with slopes in the range 1.5 < — dlogp/dlogr < 
2.5, in good agreement with the profiles observed in the nuclei of galaxies fainter than My « —20. 
However the shallow nuclear profiles of bright galaxies require a different explanation. Binary black 
holes are an inevitable result of galactic mergers, and the ejection of stars by a massive binary displaces 
a mass of order the binary's own mass, creating a core or shallow power-law cusp. This model is at 
least crudely consistent with core sizes in bright galaxies. Uncertainties remain about the effectiveness 
of stellar- and gas-dynamical processes at inducing coalescence of binary black holes, and uncoalesced 
binaries may be common in low-density nuclei. Numerical A^-body experiments are not well suited to 
probing the long-term evolution of black hole binaries due to spurious relaxation. 



1. INTRODUCTION 

The effect of a supermassive black hole on its stellar 
surroundings depends, as so often in stellar dynamics, on 
how one imagines the system evolved to its present state. 
CoUisional relaxation times are too long to have affected 
the stellar distribution in all but the densest nuclei (Faber 
et al. 1997), hence the structure and kinematics of nuclei 
are fossil relics of the interactions between stars and black 
holes. In the simplest scenario, the black hole grows by 
accreting gas on a time scale long compared with the or- 
bital periods of the surrounding stars (Peebles 1972; Young 
1980). This "adiabatic growth" model makes fairly defi- 
nite predictions about the distribution of stars near the 
black hole, predictions which are consistent with the steep 
density cusps observed in faint galaxies but which can not 
explain the flatter profiles at the centers of bright galax- 
ies. But galaxies merge, implying the formation of binary 
black holes (Begelman, Blandford & Rees 1980) which are 
efficient at displacing matter as they spiral together. The 
dynamics of black hole binaries in galactic nuclei are com- 
plex; among the unanswered questions are the long-term 
efficiency of stellar dynamical processes at extracting en- 
ergy from a binary, and whether decay of black hole bi- 
naries ever stalls at separations too great for the efficient 
emission of gravitational waves. But the binary black hole 
model is at least crudely consistent with the observed de- 
pendence of nuclear structure on galaxy luminosity. This 
article summarizes theoretical work on the single and bi- 
nary black hole models and suggests avenues for future 
progress. 

2. PRELIMINARIES 

A black hole of mass M, embedded in a galactic nucleus 
will strongly affect the motion of stars within a distance 
r = Th, the "radius of influence." A standard definition 
for rh is 
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with a the ID velocity dispersion of the stars at r 3> r/j. 
This definition had its origin in the isothermal sphere 
model for galactic nuclei; a is independent of r in such 



a model and rh is the radius at which the circular velocity 
around the black hole equals a. We now know that nuclei 
are power laws in the stellar density, p ^ , and that 
7 can lie anywhere between ~ and ~ 2.5 (Lauer, this 
volume). The velocity dispersion in a power-law nucleus 
is only constant if 7 = or 2, hence a definition like equa- 
tion is problematic. One alternative would be to define 
Th as the root of tT^(r) — GM,/r = 0. A simpler definition, 
which will be adopted in this article, is the radius at which 
the enclosed mass in stars is twice the black hole mass: 



M,{r<rh)^ 2M.. 



(2) 



This definition is exactly equivalent to equation when 
p{r) = (T^/27rGr^, the singular isothermal sphere. For an 
arbitrary power law. 
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equation ^ implies 
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Note that the "theorist's" convention is adopted here, in 
which 7 is the power-law index of the space (not projected) 
density. 

3. THE ADIABATIC GROWTH MODEL 

If a black hole grows at the center of a stellar system 
through the accretion of gas, the stellar density in the 
core will also grow as the black hole's gravity pulls in 
nearby stars (Peebles 1972; Young 1980). The change in 
the stellar density can be computed straightforwardly if it 
is assumed that the black hole grows on a time scale long 
compared with stellar orbital periods. This is reasonable, 
since even Eddington-limited accretion requires ~ 10^ yr 
to double the black hole mass, and orbital periods through- 
out the region dominated by the black hole are < 10® yr. 
Under these assumptions, the adiabatic invariants J asso- 
ciated with the stellar orbits are conserved as the black 
hole grows and the phase-space density / remains fixed 
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when expressed in terms of the J. Computing the final / 
becomes a simple matter of expressing the final orbital in- 
tegrals in terms of their initial values under the constraint 
that the adiabatic invariants remain fixed (Young 1980). 

In spherical potentials, the adiabatic invariants arc 
the angular momentum L and the radial action / — 
2 a/2 [E - $(r)] - L^/r"^, where $(r) is the gravita- 
tional potential and r± are pericenter and apocenter radii. 
It may be shown (e.g. Lynden-Bell 1963) that orbital 
shapes remain nearly unchanged when L and / are con- 
served, implying that an initially isotropic velocity distri- 
bution f{E) remains nearly isotropic after the black hole 
grows (though not exactly isotropic - see below) . The final 
/ corresponding to an initially isotropic / is then simply 

ff{Ef,L)^f,{E,,L) 

^ffiEf) (5) 

where Ef is related to Ei through the condition If{Ef, L) — 
h{E,,L). 

While the form of the nuclear density profile before the 
black hole appeared is not known, A^-body studies of struc- 
ture formation suggest that power laws are generic (Power 
et al. 2003 and references therein) and observed nuclear 
density profiles are often well described as power laws 
even on scales r :s> Vh- Setting pi oc r~''° , $i oc r'^~'^° 
(0 < 7o < 2), the initial distribution function becomes 

ME,)<^Er^, (3= (0<7o<2). (6) 

2(2 - 7o) 

To compute ff{Ef) we need a relation between Ef and 
Ei] we restrict attention to the region within the black 
hole's sphere of influence by setting Ef — t;^/2 — GM,/r. 
The radial action I{E, L) in the power-law model can not 
be computed analytically for every {E,L), but for cer- 
tain orbits Ef{Ei) has a simple form. For instance, cir- 
cular orbits have / = 0, and conservation of angular mo- 
mentum implies riMi(ri) = rfM, or r/ oc r'^~^'' ■ Thus 
Ef oc ~rj' oc -r7«-^ oc „£;(^«-4)/(2-7o)^ 

E,o^i-Ef)-^^-^"^/'^^-''°l (7) 

The same relation turns out to be precisely correct for 
radial orbits as well and is nearly correct at intermediate 
eccentricities (Gondolo & Silk 1999). Thus we can write 

ffiEf) = ME,) cx E-f' (X {-EfY, S = -^-^ (8) 

2(4 - 7o) 

and the final density profile within the sphere of infiuence 
of the black hole is 

Pf{r)= [ ff{v)(fv(x [ {^Ef^/E~~Mr) dE 

J J<S>{r) 

ocr"'^, 7 = 2 + -^. (9) 
4-70 

For < 7o < 2, 7 varies only between 2.25 and 2.5; the 
slope of the final density profile within is almost inde- 
pendent of 7o. 

The form of py(r) at r w must be computed numer- 
ically (e. g. Young 1980; Cipolhna & Berlin 1994; Cipol- 
lina 1995; Quinlan, Hernquist & Sigurdsson 1995; Gondolo 



6 Silk 1999). Figure n shows p/(r) when pi{r) oc r^^" . 
Defining r^usp to be the radius at which the inner and 
outer power laws intersect, one finds 

Tcusp^arh, 0.19<a<0.22, 0.5 < 70 < 1-5. (10) 

In early treatments of the adiabatic growth model (Pee- 
bles 1972; Young 1980), the black hole was assumed to 
grow inside of a constant-density isothermal core. The in- 
dex of the power-law cusp that forms from this initial state 
is 7 — 1.5, compared with the limiting valure 7 — 2.25 as 
7o — > in the power-law models. This difference can be 
traced to differences in the central density profile: 

p,(r) =po X (l + Cir + Car^ + ...). (11) 

The isothermal model has Ci = (an "analytic core") 
implying a phase space density that tends to a constant 
value at low energies. Other sorts of cores have Ci and 
/ diverges at low energies; for instance, the core produced 
by setting 7 = in p{r) = r~'^{l + r)~'* has f{E) 
[E — <I>(0)]^^. In fact models with finite central p's can 
be found that generate final cusp slopes anywhere in the 
range 1.5 < 7 < 2.25 (Quinlan, Hernquist & Sigurdsson 
1995). There is probably no way of ruling out an analytic 
core in the progenitor galaxy on the very small scales that 
are relevant to the later formation of a cusp, hence the 
adiabatic growth model is compatible with any final slope 
in the range 1.5 < 7 < 2.5. The upper limit could even be 
extended beyond 2.5 if 70 > 2. 

How do these predictions compare with the data? Ob- 
served luminosity profiles are well described as power laws 
at the smallest resolvable radii, and in the case of faint 
ellipticals, My > —20, the observed range of slopes is 
1-5 ^ 7 ^ 2.5 (Lauer, this volume). This is precisely the 
range in 7 predicted by the adiabatic growth model. How- 
ever in bright galaxies, 7 extends down to 0. A natural 
interpretation is that the steep cusps in faint galaxies are 
a result of adiabatic black hole growth, while some addi- 
tional mechanism, like mergers, has acted to modify the 
profiles in the brighter galaxies. 

Some fine tuning is still required to reproduce the lu- 
minosity profiles of galaxies with 1.5 < 7 < 2. These 
intermediate slopes require a shallow, core-like initial pro- 
file, and if the initial core radius exceeded ^ rh, the final 
profile will exhibit an upward inflection at r w (e.g. 
Figure 2 of Young 1980; Figure^. Such inflections are 
rarely if ever seen; observed profiles have slopes that de- 
crease smoothly inward. A way out is to require that the 
black hole mass exceeds the initial core mass so that its 
growth obliterates the core; alternatively, all nuclei with 

7 < 2 may have been the products of mergers. 

An ingenious attempt to reconcile the adiabatic growth 
model with both steep and shallow cusps was made by van 
der Marel (1999). He postulated the existence of isother- 
mal cores in the progenitor galaxies with core masses scal- 
ing as ^ L^ "'^, with L the total galaxy luminosity. Since 
M, ~ L, M,/Mcore ''^ and the black holes in faint 

galaxies would grow to dominate their cores, producing a 
cusp profile that approximates the featureless power laws 
of faint galaxies. In bright galaxies, van der Marel argued 
that the upward inflection a^t r < would be difficult 
to resolve; hence these galaxies would exhibit nearly un- 
perturbed core profiles. Van der Marel's model is intrigu- 
ing, although the assumed relation between core size and 
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Fig. 1. — Influence of the adiabatic growth of a black hole on its nuclear environment in a spherical, isotropic galaxy, (a) Density profiles 
after growth of the black hole. Initial profiles were power laws, pi oc r 'o, with 70 increasing upwards in steps of 0.25. The radial scale is 
normalized to r^i as defined in the initial galaxy (eq. EJ. The slope of the final profile at r < rj, is almost independent of the initial slope, (b) 
Velocity anisotropics after growth of the black hole. A slight bias toward circular motions appears at r < r^. 



galaxy luminosity is ad hoc, and the assumption of large, 
pre-existing cores does not fit naturally into any current 
model of galaxy formation. 

The predictions of the adiabatic growth model can be 
altered if the initial galaxy is anistropic, non-spherical or 
rotating. The simplest such case to treat is a spherical 
nucleus containing only circular orbits, with or without 
rotation (Young 1980; Quinlan, Hernquist & Sigurdsson 
1995; UUio, Zhao & Kamionkowski 2001). The derivation 
given above applies without approximations to this case: if 
the initial density profile is a power law, pi cx r~'^° , the fi- 
nal profile is also a power law with 7 = 2-1- 1/(4 — 70), 
for any 70 in the range [0,3]. As shown above, the 
isotropic model yields approximately the same final slope 
for < 7o < 2, suggesting that even extreme tangential 
anisotropies have little effect on the final profile. The ef- 
fect of radial anistropies has apparently never been inves- 
tigated; this is worth doing since real galaxies often show 
evidence for significant radial anisotropies (e.g. Kronawit- 
ter et al. 2000). 

Black hole growth in axisymmetric nuclei with analytic 
cores was considered by Leeuwin & Athanassoula (2000). 
They found little dependence of the final cusp slope on 
the degree of flattening. Merritt & Quinlan (1998) grew 
black holes of various masses in a triaxial iV-body model 
that was formed via gravitational collapse, and found p ^ 
at r < r/j. Their models evolved to axisymmetry 
after growth of the black hole, but it is now known that 
triaxiality can be maintained throughout the black hole's 
sphere of influence (Poon & Merritt 2002, 2003). Triaxial 
potentials support a wide range of different orbit families 
and the effects of black hole growth on such models have 
not been examined in detail. 

In the spherical geometry, adiabatic growth of a black 



hole induces a mild anisotropy in the stellar motions at 
r < Th due to the slightly different ways that circular 
and eccentric orbits respond to the changing potential; 
the net effect is a decrease in the average orbital eccen- 
tricity (Young 1980; Goodman & Binney 1984; Quinlan, 
Hernquist & Sigurdsson 1995; Figure ^1. If the progeni- 
tor galaxy is rotating, growth of the black hole tends to 
increase Vrot more rapidly than a (Lee & Goodman 1989; 
Leeuwin & Athanassoula 2000), although again the effect 
is slight. 

To summarize: the adiabatic growth model is limited in 
its ability to reproduce the full range of luminosity profiles 
observed in galactic nuclei. The model predicts power- 
law profiles at r < r/j with logarithmic slopes 1-5 < 7 < 
2.5. This nicely brackets the range of slopes observed in 
the nuclei of galaxies fainter than My k, —20. However 
slopes less than 7 « 1.5 are not naturally produced by the 
adiabatic growth model, and some fine tuning is required 
to avoid inflections in the profile at r « r;j. 

4. THE BINARY BLACK HOLE MODEL 

The adiabatic growth model was proposed (Peebles 
1972) before the importance of galaxy interactions and 
mergers was appreciated. We now know that supermassive 
black holes have been present in at least some spheroids 
since redshifts of 2; « 6 (e.g. Fan et al. 2001), and we 
believe that most galaxies have experienced at least one 
major merger since that time; indeed the era of peak 
quasar activity may coincide with the era of galaxy as- 
sembly via mergers (e.g. Caviliere & Vittorini 2000). If 
a nucleus forms via the merger of two galaxies contain- 
ing pre-existing black holes, the net effect on the nuclear 
density profile is roughly the opposite of what the adi- 
abatic growth model predicts: the black holes displace 
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a) 



® 



b) 







Fig. 2. — Two schemes for growing black holes at the centers of galaxies, (a) The adiabatic growth model; the stellar density within rji is 
increased as the black hole pulls in stars, (b) The binary black hole model; the inspiralling black holes displace matter within a distance 
that is roughly the separation between the two black holes when they first form a bound pair. 



matter as they spiral into the center (Figure EJ. This is 
a natural way to account for the shallow nuclear profiles 
in bright galaxies. The process may be understood as a 
sort of dynamical friction, with the "heavy particles" (the 
black holes) transferring their kinetic energy to the "light 
particles" (the stars) . However most of the energy trans- 
fer takes place after the two black holes have come within 
each other's sphere of influence, and in this regime, the 
interaction with the background is dominated by another 
mechanism, the gravitational slingshot (Saslaw, Valtonen 
& Aarseth 1974). The massive binary ejects passing stars 
at high velocities, removing them from the nucleus and 
simultaneously increasing its binding energy. 

We begin by reviewing the dynamics of massive bina- 
ries in fixed, homogeneous backgrounds, then consider the 
more difficult problem of a binary located at the center of 
an inhomogeneous and evolving galaxy. 

4.1. Dynamics of Massive Binaries 

Consider a binary system consisting of two black holes 
with masses Mi and M2, where p = M2/M1 < 1 and 
M12 = Ml + M2. Let a be the semi- major axis of the 
Keplerian orbit and e the orbital eccentricity. The binding 
energy of the binary is 



GM1M2 G/iMi2 



2a 



2a 



(12) 



with /i = M1M2/M12 the reduced mass. The relative ve- 
locity of the two black holes, assuming a circular orbit. 
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Vbin — 



GMu 



658 km s 



-1 / Ml 



VIO8M0 



1/2 



1 pc 



-1/2 



(13) 

A binary is called "hard" when its binding energy per unit 
mass, \E\/Mi2 = Gfi/2a, exceeds . For concreteness, 
a binary will here be called hard if its separation falls below 
a/i, where 



an 



Gfj_ 

4(72 



0.27pc(l +p) 



-1 / M2 



V IO^Mq 



200 km s" 

(14) 

s. Hills 1983: 



Other definitions of a/j are possible (e. 
Quinlan 1996). 

Stars passing within a distance ~ 3a of the center of 
mass of a hard binary undergo a complex interaction with 
the two black holes, followed almost always by ejection 
at velocity ~ \/ fJ,/Mi2Vbin (Saslaw, Valtonen & Aarseth 
1974). Each ejected star carries away energy and angu- 
lar momentum, causing the semi-major axis, eccentricity 
and orientation of the binary to change and the local den- 
sity of stars to drop. If the stellar distribution is assumed 
fixed far from the binary and if the contribution to the 
potential from the stars is ignored, the rate at which these 
changes occur can be computed by carrying out scatter- 
ing experiments of massless stars against a binary whose 
orbital elements remains fixed during each interaction. 

The results of the scattering experiments can be sum- 
marized via a set of dimensionless coefficients H, J, K, L, ... 
which define the mean rates of change of the parameters 
characterizing the binary and the stellar background (Hills 
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& Fullerton 1980; Roos 1981; Hills 1983, 1992; Baranov 
1984; Mikkola & Valtonen 1992; Quinlan 1996; Merritt 
2001, 2002). These coefficients are functions of the binary 
mass ratio, eccentricity and hardness but are typically in- 
dependent of a in the limit that the binary is very hard. 
The hardening rate of the binary is given by 



d /I 
dt \a 



H 



Gp 



(15) 



with p the density of stars at infinity. The mass ejection 
rate is 



dM, 



dln(l/a) 



= JMi2 



(16) 



with Afgj the mass in stars that escape the binary. The 
rate of change of the binary's orbital eccentricity is 



de 



dln(l/a) 



K. 



(17) 



The diffusion coefficient describing changes in the binary's 
orientation is 

2, r Ti* Gpa 



L- 



Mi2 a 



(18) 



with the stellar mass. Additional coefficients describe 
the rate of diffusion of the binary's center of mass, or 
"Brownian motion" (Merritt 2001). 

The binary hardening coefficient H reaches a constant 
value of ~ 16 in the limit a ^ a^, with a weak dependence 
on M2/M1 (Hills 1983; Mikkola & Valtonen 1992; Quin- 
lan 1996). In a fixed background, equation (|15|) therefore 
implies that a hard binary hardens at a constant rate: 



a[t) ah a 



tfi) , t> th, a{th) = ah- (19) 



If the supply of stars remains steady, hardening continues 
at this rate until the components of the binary come close 
enough together that the emission of gravitational radia- 
tion is important. In this regime, gravity wave coalescence 
takes place in a time: 



F(e) = (l-e2) 



-7/2 



73 
— f 

24 



37 4 
— e 
96 



(20) 



(Peters 1964). Coalescence in a time tgr occurs when a — 



Qgr, where 



P' 



,3/4 



j,gr 



3/2 V 200 km s" 



-7/8 



/ tgr \ 



-1/4 



Vlov / 

(21) 

The M, — a relation has been used to express M\i in terms 
of a. For mass ratios p of order unity and e « 0, the binary 
must decay by a factor of ~ 10^ in order for gravitational 
radiation to induce coalescence in a time shorter than 10^ 
yr. Less decay is required if the binary is eccentric or if 
M2 < Ml. 

If the binary manages to shrink by such a large factor, 
the damage done to its stellar surroundings will be consid- 
erable. The mass ejected by the binary in decaying from 
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Fig. 3. — Mass ejected by a decaying binary, in units of M12 = 
Ml -A- Ml (solid lines) or Mi (dashed lines). Curves show mass that 
must be ejected in order for the binary to reach a separation where 
the emission of gravity waves causes coalescence on a time scale of 
10^" yr (lower), lO'' yr (middle) and 10* yr (upper). 



o/i to is given by the integral of equation (|16|l : 



A/12 



J{a) 



da. 



(22) 



Figure 1.2 shows M^j as a function of the mass ratio 
M2/M1 for a = 200 km s~^ and various values of tgr- 
The mass ejected in reaching coalescence is of order Mi 2 
for equal-mass binaries, and several times M2 when M2 <C 
Ml. A black hole that grew to its current size through 
a succession of mergers should therefore have displaced a 
few times its own mass in stars. If this mass came mostly 
from stars that were originally in the nucleus, the den- 
sity within r/j would drop drastically and the hardening 
would slow. Without some way of replenishing the supply 
of stars, decay could stall at a separation much greater 
than agr- This is the "final parsec problem": how to avoid 
stalling and bring the black holes from their separation of 
^ 1 pc, when they ffist form a hard binary, to ~ 10"^ pc, 
where gravity wave emission is efficient. 

Changes in the binary's orbital eccentricity (equation 
I17f) are potentially important because the gravity wave 
coalescence time drops rapidly as e ^ 1 feauation l20l) . For 
a hard binary, scattering experiments give K{e) ~ Koe{l~ 
e^), with Ko « 0.5 for an equal- mass binary (Mikkola & 
Valtonen 1992; Quinlan 1996). The dependence of K on 
M2/M1 is not well understood. The implied changes in e 
as a binary decays from a = to agr are modest, Ae < 
0.2, for all initial eccentricities. The rms change in the 
orientation of the binary's spin axis feauation ll8|) is SO « 
2(m,/Afi2)^/2 log^/2(a/i/a), which is of order one degree or 
less if m, « Mq. Reorientations of the binary's angular 
momentum affect the spin direction of the coalesced black 
hole and the direction of any associated jet (Merritt 2002). 

4.2. Binary Black Holes in Galaxies 

The scattering experiments summarized above treat the 
binary's environment as fixed and homogeneous. In real- 
ity, the binary is embedded at the center of an inhomo- 
geneous and evolving galaxy, and the supply of stars that 
can interact with it is limited. 
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In a fixed spherical galaxy, stars can interact with the 
binary only if their pericenters lie within ~ 7?, x a, where 
7^ is of order unity. Let Lie = 'R-a^y2[E - <^{TZa)] « 
\/2GMi2R-a^ the angular momentum of a star with peri- 
center Tla. The "loss cone" is the region in phase space 
defined by L < Lie- The mass of stars in the loss cone is 



dE / dL^4Tr\f{E,L'')P{E,L^) 



« 87^2GMl2TO*7^a J dEf{E)Prad{E). (23) 

Here P is the orbital period; in the final line, / is assumed 
isotropic and P has been approximated by the period of a 
radial orbit of energy E. An upper limit to the mass that 
is available to interact with the binary is ~ Midah), the 
mass within the loss cone when the binary first becomes 
hard; this is an upper limit since some stars that are ini- 
tially within the loss cone will "fall out" as the binary 
shrinks. Assuming a singular isothermal sphere for the 
stellar distribution, p cx r"^, and taking the lower limit of 
the energy integral to be ^(ah), equation (|23|l implies 



Miciah) ~ STZp.. 



(24) 



We can compute the change in a that would result if the 
binary interacted with this entire mass, by using the fact 
the mean energy change of a star interacting with a hard 
binary is ~ 3GyLt/2a (Quinlan 1996). Equating the en- 
ergy carried away by stars with the change in the binary's 
binding energy gives 



3Gp 



dM 



or 



In 



(t) 



, AM 9TZn 



Mr 



Ml 



GMiM2 ^n 
2 la 



97^ 



(25) 



p EE Ma /Ml 



(26) 

if AM is equated with Mic- Only for very low mass ratios 
(M2 < 10~^Mi) is this decay factor large enough to sat- 
isfy equation (|21l) . but the time required for such a small 
black hole to reach the nucleus is likely to exceed a Hubble 
time (Merritt 2000). Hence even under the most favorable 
assumptions, the binary would not be able to interact with 
enough mass to reach gravity- wave coalescence. 

But the situation is even worse than this, since not all of 
the mass in the loss cone will find its way into the binary. 
The time scale for the binary to shrink is comparable with 
stellar orbital periods, and some of the stars with rp^ri ~ 
ah will only reach the binary after a has fallen below ~ ah- 
We can account for the changing size of the loss cone by 
writing 



1 dM, 



dM 

.lE,(t)P{E) dE 



dE 



= 8TT^GMi2rru'Jla{t) / ME)dE, 

JEoit) 



(27) 



where M{t) is the mass in stars interacting with the bi- 
nary and fi(E) is the initial distribution function; setting 
P{Eq) = t reflects the fact that stars on orbits with pe- 
riods less than t have already interacted with the binary 
and been ejected. Combining equations H25|l and (|27|l . 



d_ 

di 



- ) « 2ATT^'R.Gm^ [ ME)dE. (28) 



Solutions to equation (|28|) show that a binary in a sin- 
gular isothermal sphere galaxy stalls at ah/ a « 2.5 for 
M2 = Ml, compared with ah/ a ~ 10 if the full loss cone 
were depleted (equation I26|l . In galaxies with shallower 
central cusps, decay of the binary would stall at even 
greater separations. 

It is therefore entirely possible that uncoalesced binaries 
exist at the centers of many galaxies, particularly galaxies 
like large ellipticals with low central densities. But there 
is circumstantial evidence that long-lived black hole bina- 
ries are rare. Some radio galaxies show clear evidence of 
a recent flip in the black hole's spin orientation, as would 
occur when two black holes coalesce (Dennett-Thorpe et 
al. 2002), and their numbers are consistent with a coa- 
lescence rate that is roughly equal to the galaxy merger 
rate (Merritt & Ekers 2002). The almost complete lack of 
correlation between jet directions in Seyfert galaxies and 
the angular momenta of their disks (Ulvestad & Wilson 
1984; Kinney et al. 2000) also suggests that black hole 
coalescences were common in the past. There are few if 
any "smoking gun" detections of binary black holes among 
AGN (e.g. Halpern & Eracleous 2000); the best, but stiU 
controversial, case is OJ 287 (Pursimo et al. 2000). 

Below are discussed some additional mechanisms that 
have been proposed for extracting energy from binary 
black holes and hastening their decay. When the agent in- 
teracting with the binary is stars, continued decay of the 
binary implies continued destruction of the stellar cusp. 
However for most of the mechanisms discussed below, the 
detailed effects on the stellar distribution have yet to be 
worked out. 

4.2.1. Scattering of stars into the loss cone 

Destruction of a pre-existing stellar cusp generates strong 
gradients in the phase-space density at L Lie, the an- 
gular momentum of an orbit that lies just outside the loss 
cone. A small perturbation can deflect a star on such 
an orbit into the loss cone. This process has been stud- 
ied in detail in the context of gravitational scattering of 
stars into the tidal disruption sphere of a single black hole 
(Lightman & Shapiro 1977; Cohn & Kulsrud 1978). Once 
a steady-state flow of stars into the loss cone has been set 
up, the distribution function near Lie has the form: 



(29) 



(Lightman & Shapiro 1977), where R = L^/Ll{E), Lc{E) 
is the angular momentum of a circular orbit of energy E, 
and / is the distribution function far from the loss cone, 
assumed to be isotropic. The mass flow into the central 
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object is TO* / T{E)dE, where 
HE)dE^,.^LliE)[j'^J^^ 



2R 



f 



ln(l/i?;c) 



-dE. 



^^^^ 

The quantity in brackets is the orbit-averaged diffusion 
coefficient in R. Yu (2002) evaluated the contribution of 
two-body scattering to the decay rate of a massive binary 
and found that it was usuaUy too smaU to overcome the 
stalhng that occurs when the loss cone is first emptied. 

Standard loss cone theory as applied by Yu (2002) as- 
sumes a quasi-steady-state distribution of stars in phase 
space near Lie. This assumption is appropriate at the cen- 
ter of a globular cluster, where relaxation times are much 
shorter than the age of the universe, but is less appropri- 
ate for a galactic nucleus, where relaxation times almost 
always greatly exceed a Hubble time (e.g. Faber et al. 
1997). The distribution function f{E,L) immediately fol- 
lowing the formation of a hard binary is approximately a 
step function. 



(31) 



much steeper than the ^ In i dependence in a collisonally 
relaxed nucleus (equation 1291) . Since the transport rate 
in phase space is proportional to the gradient of / with 
respect to L, steep gradients imply an enhanced flux into 
the loss cone. The total mass in stars consumed by the 
binary can exceed the predictions of the standard model 
by factors of a few, implying greater cusp destruction and 
more rapid decay of the binary (Milosavljevic & Merritt 
2003b). This time-dependent loss cone refilling might be 
particularly effective in a nucleus that that continues to 
experience mergers or accretion events, in such a way that 
the loss cone repeatedly returns to an unrelaxed state with 
its associated steep gradients. 

4.2.2. Re-ejection 

Unlike the case of tidal disruption of stars by a single 
black hole, a star that interacts with a massive binary re- 
mains inside the galaxy and is available for further interac- 
tions. In principle, a single star can interact many times 
with the binary before being ejected from the galaxy or 
falling outside the loss cone; each interaction takes addi- 
tional energy from the binary and hastens its decay. Con- 
sider a simple model in which a group of N stars in a 
spherical galaxy interact with the binary and receive a 
mean energy increment of {^E). Let the original energy 
of the stars be Eq. Averaged over a single orbital period 
P{E), the binary hardens at a rate 



'GA/iA/2\ _ N{AE) 
dt V 2a J P{E) 



(32) 



In subsequent interactions, the number of stars that re- 
main inside the loss cone scales as Lf^ cx a while the ejec- 
tion energy scales as ~ a""'^. Hence N{AE) cx a^a~^ cx a°. 
Assuming the singular isothermal sphere potential for the 
galaxy, one finds 



ah 
a{t) 



1 



In 



1 



m.^N(AE) t-th 



2na^ P{Eo) 



(33) 



(Milosavljevic & Merritt 2003b). Hence the binary's bind- 
ing energy increases as the logarithm of the time, even 
after all the stars in the loss cone have interacted at least 
once with the binary. Re-ejection would occur differently 
in nonspherical galaxies where angular momentum is not 
conserved and ejected stars could miss the binary on their 
second passage. However there will generally exist a sub- 
set of orbits defined by a maximum pericenter distance 
< a and stars scattered onto such orbits can continue to 
interact with the binary. 

4.2.3. Chaotic loss cones 

Loss cone dynamics are qualitatively different in non- 
axisymmetric (triaxial or bar-like) potentials, since a much 
greater number of stars may be on "centrophilic" orbits 
which take them near to the black hole(s) (Gerhard & 
Binney 1985). Triaxial models need to be taken seriously 
given recent demonstrations (Poon & Merritt 2002, 2003) 
that black hole nuclei can be stably triaxial and that the 
fraction of mass on centrophilic - typically chaotic - orbits 
can be large. The frequency of pericenter passages, rpori < 
d, for a chaotic orbit of energy Erne, triaxial black-hole 
nuclus is roughly linear in d, N{rpci-i < d) ~ d x A{E) 
(Merritt & Poon 2003). The total rate at which stars pass 
within a distance TZa of the massive binary is therefore 



dM 
IT 



TZa / A{E)Mc{E)dE 



(34) 



where Mc{E) is the mass on centrophilic orbits. The im- 
plied feeding rate can be comparable to that in a spherical 
nucleus with a constantly-refilled loss cone, and orders of 
magnitude greater than in a loss cone that is re-supplied 
via star-star interactions (Merritt & Poon 2003). Even 
transient departures from axisymmetry, for instance dur- 
ing mergers, might result in substantial loss cone refilling 
due to this mechanism. 

4.2.4. Multiple black holes 

If binary decay stalls, an uncoalesced binary may be 
present in a nucleus when a third black hole, or a second 
binary, is deposited there following a subsequent merger. 
The multiple black hole system that forms will engage 
in its own gravitational slingshot interactions, eventually 
ejecting one or more of the black holes from the nucleus 
(though probably not from the galaxy). This process has 
been extensively modelled assuming a fixed potential for 
the galaxy (e.g. Valtaoja, Valtonen & Byrd 1989; Mikkola 
& Vahonen 1990; Vaffonen ct al. 1994). The effect on 
the stellar distribution of iV > 2, interacting black holes 
is not well understood, though iV-body simulations with 
10 — 20 massive particles and a "live" background show 
that the black holes displace ~ 10 times their own mass 
before being ejected (Merritt & Milosavljevic 2003). The 
separation and eccentricity of the dominant binary can 
change dramatically during each interaction and this may 
be an effective way to shorten the gravity wave coalescence 
time. In a wide, hierarchical triple, the eccentricity of the 
dominant binary oscillates through a maximum value of 



5cos2i/3, |cosi| < \/3/5, with i the mutual in- 
clination angle (Kozai 1962). Blaes, Lee & Socrates (2002) 
estimate that the coalescence times in equal-mass, hierar- 
chical triples can be reduced by factors of ^ 10 via the 
Kozai mechanism. 
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4.2.5. Gas dynamics 

If the inner '--^ 1 pc of the nucleus contains a mass in 
gas comparable to M2, torques from the gas will cause the 
orbit of the smaller black hole to decay in a time of or- 
der the gas accretion time (Syer & Clarke 1995; Ivanov, 
Papaloizou & Polnarev 1999). Given standard assump- 
tions about accretion disk viscosities, the gas-dynamical 
decay rate would exceed that from gravity wave emission 
for a > aacc, where 

— - 1 X 10"' drrf' ( (35) 

ah VO.l/ V200 kms"V 

(Armitage & Natarajan 2002). If the orbit of the sec- 
ondary is strongly inclined with respect to the accretion 
disk around the larger black hole, its passages through the 
disk could generate periodic outbursts, and this has been 
suggested as a model for the ~ 12 yr cycle of optical flaring 
observed in the blazar OJ 287 (Lehto & Valtonen 1996). 
Gas deposition of the required magnitude almost certainly 
occurred during the quasar epoch, although it is less clear 
that this mechanism is effective for galaxies in the current 
universe. 

5. TV-BODY STUDIES 

Unless great care is taken, A^-body studies of binary 
black hole dynamics are unlikely to give an accurate pic- 
ture of the evolution expected in real galaxies. This follows 
from the result (§ 14.21) that time scales for two-body scat- 
tering of stars into the binary's loss cone are of order the 
Hubble time or somewhat longer in real galaxies. In N- 
body simulations, relaxation times are shorter by factors 
of ^ A^/10^^ than in real galaxies, hence the long-term evo- 
lution of the binary is likely to be dominated by spurious 
loss cone refilling, Brownian motion of the black holes and 
other noise-driven effects (Milosavljevic & Merritt 2003b). 
The stalling that is predicted in the absence of loss cone 
refilling in real galaxies (Figure 0J can only be reproduced 
via iV-body codes if the mean field is artificially smoothed 
and the black hole binary is "nailed down" (e.g. Quinlan 
& Hernquist 1997). 

A^-body studies are most useful at characterizing the 
early stages of binary formation and decay, or simulating 
the disruptive effects of a single black hole on an infalling 
galaxy; indeed scattering experiments (§ I4.1|l are almost 
useless in these regimes. Due to algorithmic limitations, 
most A'-body studies (e.g. Ebisuzaki, Makino & Okumura 
1991; Makino et al. 1993; Governato, Colpi & Maraschi 
1994; Makino & Ebisuzaki 1996; Makino 1997; Nakano 
& Makino 1999a, b; Hemsendorf, Sigurdsson & Spurzem 
2002) have been based on galaxy models with unrealis- 
tically large cores. A^-body merger simulations using re- 
alistically dense initial conditions (HoUey-Bockclmann & 
Richstone 2000; Merritt & Cruz 2001; Merritt et al. 2002) 
show that the black hole in the larger galaxy is efficient 
at tidally disrupting the steep cusp in the infalling galaxy, 
producing a remnant with only slightly higher central den- 
sity than that of the giant galaxy initially. This result 
helps explain the absence of dense cusps in bright galax- 
ies (Forbes, Franx & lUingworth 1995) and the persistence 
of the "core fundamental plane" in the face of mergers 
(HoUey-Bockelmann & Richstone 1999). 



Quinlan & Hernquist (1997) studied the evolution of a 
black hole binary inside cuspy models with p ^ r^^ and 
p ^ and a range of black hole masses and particle 
numbers, N < 2 x 10^. Their A^-body code was unable to 
simulate an actual merger and all of their detailed results 
were derived from initial conditions consisting of a sin- 
gle galaxy into which two "naked" black holes were sym- 
metrically dropped. This artificial starting configuration 
produced substantial evolution of the cusp before the for- 
mation of the binary. The late evolution of the binary 
was found to be strongly dependent on A^, due in part 
to spurious Brownian motion of the black hole particles. 
The cores that formed were characterized by strong veloc- 
ity anisotropics, at 3> cr^, due to the ejection of stars on 
eccentric orbits. 

Milosavljevic & Merritt (2001) followed the evolution 
of cuspy {p ^ r~^) galaxy models containing black holes, 
starting from pre-merger initial conditions and continuing 
until the binary separation had decayed a factor of 10 
below ah- The initially steep nuclear cusps were converted 
to shallower, p ~ profiles shortly after the black holes 
had formed a hard binary; thereafter the nuclear profile 
evolved slowly toward shallower slopes. The decay rate of 
the binary was found not to be strongly dependent on A^, 
probably due to the fact that the loss cone was contin- 
uously refilled by two-body scattering in their coUisional 
simulation (Milosavljevic & Merritt 2003b). The veloc- 
ity anisotropics created during formation of the core were 
much milder than in the simulations of Quinlan & Hern- 
quist (1997), similar in magnitude to the anisotropics pre- 
dicted by the adiabatic growth model (Figure^. 

An A^-body test of the time-dependent loss cone dynam- 
ics described in §4.2 would be possible using a coUisional 
code with A^ > 10^ (Milosavljevic & Merritt 2003b). Such 
large particle numbers demand a novel A^-body algorithm 
coupled with special-purpose hardware like the GRAPE-6. 

6. OBSERVATIONAL EVIDENCE FOR THE BINARY BLACK 
HOLE MODEL 

Figure|21shows that a massive binary must eject of order 
its own mass in reaching a separation of agr if M2 ~ Mi , 
or several times M2 if M2 <C Mi. These numbers should 
be interpreted with caution since: (1) binaries might not 
decay as far as agr, or the final stages of decay might be 
driven by some process other than energy exchange with 
stars; (2) the definition of "ejection" used in Figure 13 is 
escape of a star from an isolated binary, and does not 
take into account the confining effect of the nuclear po- 
tential; (3) the effect of repeated mergers on nuclear den- 
sity profiles, particularly mergers involving very unequal- 
mass binaries, is poorly understood. Nevertheless, a clear 
prediction of the binary black hole model is that galactic 
mergers should result in the removal of a mass of order 
Afi2 from the nucleus. As in the adiabatic growth model, 
we are handicapped in testing the theory by lack of knowl- 
edge of the primordial nuclear profiles. A reasonable guess 
is that all galaxies originally had steep, p ~ density 
cusps, since these are generic in the low-luminosity ellip- 
ticals which are least likely to have been influenced by 
mergers, and since the adiabatic growth model predicts 
p (Figure P). 

The "mass deficit" is defined as the difference in inte- 
grated mass between the observed density profile, and a 
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M. M. M. 

Fig. 4. — Mass deficit versus black hole mass for three different assumed values of 70, the logarithmic slope of the density cusp before 
energy input from the black holes, (a) 70 = 2; (b) 70 = 1.75; (c) 70 = 1.5. Units arc solar masses. Solid lines are M^^j = M, (Milosavljevic 
et al. 2002). 



p cc r ^° profile extrapolated inward from the break radius 



Mdef = 47r 



Pin) - 
n 



-70 



p(r) 



r^dr. (36) 



Milosavljevic et al. (2002) and Ravindranath, Ho & Filip- 
penko (2002) computed Mdef in samples of "core" -profile 
elliptical galaxies. The former authors found a good pro- 
portionality between M^e/ and 7\f,, with {Mdef /M,) ~ 1 
for 70 = 1.5 and (Mdef/M,) w 10 for 70 = 2 (Fig- 
ure01). These numbers are within the range predicted by 
the binary black hole model, given the uncertainties asso- 
ciated with the effects of multiple mergers. Ravindranath, 
Ho & Filippenko (2001) computed black hole masses us- 
ing a shallower assumed slope for the M, — a relation, 
M, oc a^'"^^, and found a steeper, nonlinear dependence of 
Mdef on M,. 

More rigorous tests of the binary black hole model will 
require a better understanding of the expected form of 
p{r). As discussed above, while the best current iV-body 
simulations suggest p ~ following binary formation 
(Milosavljevic & Merritt 2001), the simulations are domi- 
nated by noise over the long term. If the decay stalls, the 
predicted density profile can be very different: a "hole" 
forms inside of ^ Sastaii (e.g. Figure 1 of Zier & Bier- 
mann 2001). Central minima may in fact have been seen 
in the luminosity profiles of a few galaxies (Lauer et al. 
2002). 

Other processes could result in energy exchange between 
binary black holes and stars. If the binary eventually 
coalesces, the gravitational radiation carries a linear mo- 
mentum leading to a recoil of the coalesced black hole at 
a velocity VrecoH ^ 10^ — 10^^ km (Bekenstein 1973; 
Fitchett 1983; Eardley 1983), and possibly even higher if 
the black holes were rapidly spinning prior to coalescence 
(Redmount & Rees 1989). Recoil velocities of this order 
would eject the black hole from the center of the nucleus 
and its subsequent infall would displace stars. Quanti- 
tative evaluation of this effect will require more accurate 
estimates of frecoii based on fully general-relativistic cal- 



culations of black hole coalescence. 

Equating Mdef with the mass removed by the binary 
black hole is only justified if galaxy mergers leave nuclear 
density profiles nearly unchanged in the absence of the 
gravitational slingshot. This is known to be the case in 
equal-mass mergers between galaxies with power-law cusps 
(Barnes 1999; Milosavljevic & Merritt 2001), though less 
is known about the effects of unequal-mass mergers. It is 
difficult to see how the break radii in galaxies with ^ rh 
(e.g. NGC 3640, 4168) can be accounted for by the binary 
black hole model unless mergers sometimes induce changes 
in p[r) well beyond the radius of influence of the black 
holes. 

Another unsolved problem is the persistence of steep, 
p ^ r^"^ nuclear profiles in fainter galaxies; even these 
galaxies should have undergone at least one merger since 
their formation. One possibility is that faint ellipticals ex- 
perienced their last major merger at a time when spheroids 
contained a much larger gas fraction than at present (e.g. 
Kauffmann & Haehnelt 2000), and that binary decay was 
driven by gas dynamics during this last merger. 

A major focus of future work should be to calculate the 
evolution of p{r) as predicted by the various scenarios for 
binary decay discussed in this article. 



An abridged version of this article will appear in Carnegie 
Observatories Astrophysics Series, Vol. 1: "Coevolution 
of Black Holes and Galaxies," edited by L. C. Ho (Cam- 
bridge: Cambridge Univ. Press). I thank M. Milosavlje- 
vic for his detailed reading of the manuscript. This work 
was supported by NSF grants AST 00-71099 and AST 02- 
06031, by NASA grants NAG5-6037 and NAG5-9046, and 
by grant HST-AR-08759 from STScI. 
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